Bilayer graphene: gap tunability and edge properties 



Eduardo V Castro^ N M R Peres^, J M B Lopes dos Santos^ 
F Guinea^ and A H Castro Neto^ 

^ CFP and Departamento de Fi'sica, Faculdade de Cincias Universidade do Porto, P-4169-007 
Porto, Portugal 

^ Center of Physics and Departamento de Fi'sica, Universidade do Minho, P-4710-057 Braga, 
Portugal 

^ Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco, E-28049 Madrid, Spain 
* Department of Physics, Boston University, 590 Commonwealth Avenue, Boston, MA 
02215,USA 

E-mail: evcastro@fc.up.pt 

Abstract. Bilayer graphene - two coupled single graphene layers stacked as in graphite - 
provides the only known semiconductor with a gap that can be tuned externally through electric 
field effect. Here we use a tight binding approach to study how the gap changes with the applied 
electric field. Within a parallel plate capacitor model and taking into account screening of the 
external field, we describe real back gated and/or chemically doped bilayer devices. We show 
that a gap between zero and midinfrared energies can be induced and externally tuned in these 
devices, making bilayer graphene very appealing from the point of view of applications. However, 
applications to nanotechnology require careful treatment of the effect of sample boundaries. This 
being particularly true in graphene, where the presence of edge states at zero energy - the Fermi 
level of the undoped system - has been extensively reported. Here we show that also bilayer 
graphene supports surface states localized at zigzag edges. The presence of two layers, however, 
allows for a new type of edge state which shows an enhanced penetration into the bulk and 
gives rise to band crossing phenomenon inside the gap of the biased bilayer system. 



1. Introduction 

The recent production of graphene [1], the first truly one-atom thick material, enabled several 
unthoughtful experiments where charge carriers were undoubtedly shown to be massless with a 
linear dispersion relation [2]. 

In addition to single layer graphene (SLG), few-layer graphene can also be isolated. Of 
particular interest to us is the double layer graphene system, where two carbon layers are placed 
on top of each other according to the usual Bernal Ai?-stacking. The low-energy properties of 
this so-called bilayer graphene (BLG) are then described by massive Dirac fermions [2], with a 
quadratic dispersion close to the neutrality point and a Dirac fermion mass originating from the 
inter-plane hopping energy t^. 

One of the most remarkable properties of BLG is the ability to open a gap in the spectrum 
by electric field effect - biased BLG. This has been realized experimentally, providing the 
first semiconductor whose gap can be tuned externally O [U [5]. In the absence of external 
perpendicular electric field - unbiased BLG - the system is characterized by four bands, two 
of them touching each other at zero energy, and giving rise to the massive Dirac fermions 
mentioned above, and other two separated by an energy ibt_L. Hence, an unbiased BLG is 



a two-dimensional zero-gap semiconductor. At the neutrality point the conductivity shows a 
minimum of the order of the conductance quantum [6], a property shared with SLG [2]. This 
prevents standard logic applications where the presence of a finite gap producing high on-off 
current ratios is of paramount importance. The fact that a simple perpendicular electric field is 
enough to open a gap, and even more remarkable, to control its value, clearly demonstrates the 
potential of this system for carbon-based electronics. A new plateau at zero Hall conductivity 
[3] and ferromagnetic behavior at low densities are some of the gap consequences [7j. Moreover, 
the induced gap has been shown to be robust to the presence of disorder |8, S\ ■ 

Although such a tunable gap by the electric field effect is not possible in SLG, band gaps 
can still be engineered by confining graphene electrons in narrow ribbons [1U\ However, the 
lateral confinement brings about the presence of edges, which in graphene can have profound 
consequences on electronics. This is essentially due to the rather different behavior of the 
two possible (perfect) terminations in graphene: zigzag and armchair. While zigzag edges 
support localized states, armchair edges do not [12j. These surface (or edge) states occur at zero 
energy, the same as the Fermi level of undoped graphene, meaning that low energy properties 
may be substantially altered by their presence. The self-doping phenomenon |13j and the edge 
magnetization with consequent gap opening jl2] are among edge states driven effects. In BLG, 
the question regarding the presence of edge states is also pertinent. Firstly, zigzag edges are 
among the possible terminations in BLG, and secondly, the presence of edges is unavoidable in 
tiny devices. 

The paper is organized as follows: in Sec.[2]the lattice structure of BLG and the tight-binding 
Hamiltonian are presented; bulk electronic properties are discussed in SecO the presence of edge 
states localized at zigzag edges of BLG is discussed in Sec. SJ Sec. [5] contains our conclusions. 

2. Model 

The lattice structure of BLG is shown in Fig. [^a). Here we consider only AB-Jierna\ stacking, 
where the top layer has its A sublattice on top of sublattice B of the bottom layer. We use 
indices 1 and 2 to label the top and bottom layer, respectively. The basis vectors may be written 
as ai = aex and a2 = a(e^ — y/3ey)/2, where a = 2.46A. 

In the tight-binding approximation, the in-plane hopping energy, t ~ 3 eV, and the inter-layer 
hopping energy, t±/t ~ 0.1 <^ 1 [2], define the most relevant energy scales (see Fig. [1]). The 
simplest tight-binging Hamiltonian describing non-interacting vr— electrons in BLG then reads: 
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H = -^Y^Yl [<a(r)^*,<x(r) + aljr)bi,^{r - ai) + aljr)bi,^{r - aa) + h.c] 

r,(T i=l 

-t±_ ^ [4^^(r)62,a(r) + h.c] +-^^ [nAiir) + nBi(r) - nA2(r) - nB2(r)] , (1) 

r,cr r,cr 

where aj^(j(r) [6i.CT(r)] is the annihilation operator for electrons at position r in sublattice Ai 
(Bi), i = 1,2, and spin a, and ra^j(r) and nsi(r) are number operators. We are interested in 
the properties of BLG in the presence of a perpendicular electric field - the biased BLG. The 
external perpendicular electric field gives rise to an electrostatic energy difference between the 
two layers, which we parametrize by V per electron. The effect of this energy difference between 
layers is accounted for by the last term in Eq. ([1]). 

3. Bulk electronic properties 

In the absence of applied perpendicular electric field the system has four bands given by 



E^^ = ±^Jel + tl/4±tj2, 



(2) 




Figure 1. (a) Lattice structure of a graphene bilayer. (b)-(c) Band structure near K (or K') 
for y = and V = t±, respectively. We used t± = 0.2t. (d) Variation of the gap Ag with V for 
t± = {0.1,0.2,0.3}t. 

where is the dispersion of SLG, = + 2cos{akx) + 4 cos {akx/2) cos (a/cy \/3/2)] . The 
band structure defined by Eq. ([2]) is shown in Fig. mb). In undoped BLG the system has exactly 
one electron per vr orbital, that is, it is a half-filled system. In this case the chemical potential 
crosses exactly at the K and K' points (the Dirac points) at the corners of the first Brillouin 
zone. As can be seen in Fig. [T](b), at these corners the dispersion is parabolic, E(j)) ~ ±Vpp'^ /t±, 
with p = hq, where q measures the distance in momentum space relatively to the Dirac points, 
and vp = \piathr^ l2 is the SLG Fermi velocity. Thus, low-energy quasiparticles in BLG are 
massive, with a light effective mass given by m* = i±/(2up) ~ O.OSme, where me is the bare 
electron mass [2]. Although this behavior is clearly different from the low-energy massless 
Dirac fermions found in SLG, e(p) ~ ivpp, it has been shown that the effective 2-band model 
describing low-energy physics in BLG is still a pseudo-spin Hamiltonian [6]. Consequently, BLG 
is also far from a standard two-dimensional electron gas; the observed anomalous quantum Hall 
effect is an example of such an unconventional behavior |6j . 

Now we address the electronic structure of the biased BLG. The spectrum of Eq. ([1]) for 
y / reads: 

Ei^{V)=±s^el + 4/2 + ^2/4 ± + {t\ + V^)\. (3) 

The resulting band structure is shown in Fig. mc). As can be seen from Eq. the V = ^ 
gapless system turns into a semiconductor whose gap is controlled by V . The gap between 
conduction and valence bands, A^, is given by 

A, = ^t\vy{t\^vy (4) 

Note that V parametrizes the effect of a perpendicular electric field, and therefore can be 
controlled externally. This means, as a consequence of Eq. (jH, that the biased BLG provides 
a semiconductor with a gap that can be tuned externally by electric field effect. From Eq. @ 
it can be seen that for both V -^t^ and V ^ t one finds A^ ~ V . However, there is a region 
foi" ^± ^ ^ ^ 6t where the gap shows a plateau Ag ~ as depicted in Fig. HKd). The plateau 
ends when 1/ ~ 6i (not shown). 

So far we have considered y as a band parameter that controls the gap. However, the 
parameter V can be related with the perpendicular electric field applied to BLG, avoiding the 
introduction of an extra free parameter in the present theory. If E = Ee^ is the perpendicular 
electric field felt by electrons in BLG, the corresponding electrostatic energy U{^z) for an electron 
of charge — e is related to the electric field as eE = dU{z)/dz, and thus V is given by 

V = U{zi) - U{z2) = eEd, (5) 

where zi and Z2 is the position of layer 1 and 2, respectively, and d = zi — Z2 = 3.4 A is the 
inter-layer distance. Given the experimental conditions, the value of E can be calculated under 
a few assumptions, as detailed in the following. 



3.1. Real biased bilayer devices 

If we assume the electric field E in Eq. ([5|) to be due exclusively to the external electric field 
applied to BLG, E = Egxt, all we need in order to know V is the value of Egxt- V = eE^xtd. The 
experimental realization of a biased BLG has been achieved in epitaxial BLG through chemical 
doping [3] and in back gated exfoliated BLG OS]. In either case the value of E^xt can be 
extracted assuming a simple parallel plate capacitor model. 

In the case of exfoliated BLG, devices are prepared by micromechanical cleavage of graphite 
on top of an oxidized silicon wafer (t = 300 nm of Si02)- A back gate voltage Vg applied between 
the sample and the Si wafer induces charge carriers due to the electric field effect, resulting in 
carrier densities ng = aVg relatively to half-filling. The geometry of the resulting capacitor 
determines the coefficient a = esi02^o/(ei) — 7.2 x 10^'^cm~^/V, where £3102 = 3.9 and eo are 
the permittivities of Si02 and free space, respectively. In order to control independently the 
gap value and the Fermi level, in Ref. [5] the devices have been chemically doped by deposition 
of ammonia (NH3) on top of the upper layer, which adsorbed on graphene and effectively acted 
as a top gate providing a fixed electron density tiq. Charge conservation then implies a total 
density n in BLG given hy n = Ug + UQ. Extending the parallel plate capacitor model to include 
the effect of dopants, the external field E^xt is the result of charged surfaces placed above and 
below BLG. The accumulation or depletion layer in the Si wafer contributes with an electric 
field E'b = eng/(2ereo)) while dopants above BLG effectively provide the second charged surface 
with electric field Et = —enQ/{2er£o), where Sr is the bilayer relative dielectric constant. Adding 
the two contributions, Egxt = Ey, + Et, and making use of the charge conservation relation, we 
arrive at an electrostatic energy difference V that depends linearly on the BLG density. 



In the case of epitaxial BLG, devices are grown on SiC by thermal decomposition. Due to 
charge transfer from SiC substrate to film, the as-prepared BLG devices appear electron doped 
with density Tia. The substrate's depletion layer provides the external electric field necessary to 
make the system a biased BLG. In Ref. [3] the BLG density n was varied by doping the system 
with potassium (K) on top of the upper layer, which originates an additional charged layer 
contributing to the external electric field. Applying the same parallel plate capacitor model as 
before, we get an electrostatic energy difference that can be written as 



In the inset of Fig. ^a) we compare Eq. ([7]) with experimental results for V from Ref. [3j, 
obtained by fitting angle resolved photo-emission spectroscopy (ARPES) measurements. For this 
particular biased BLG realization, the as-prepared carrier density was na ~ 10 x 10^^ cm^^. From 
Eq. ([7]), this ria value implies a zero V for the bilayer density n^^ ~ 20 x 10^^ cm~^, and therefore 
zero gap through Eq. Experimentally, a zero gap was found around n^^^ ~ 23 x 10^^ cm~^. 
Given the simplicity of the theory, it can be said that n*^ and n'^^^ are in good agreement. 
However, the agreement is only good at V = 0, since the measured V is not a linear function of 
n, as Eq. ([7|) implies. This is an indication that screening due to correlation effects should be 
taken into account. 

3.2. Screening correction and gap behavior 

In deriving Eqs. ([6]) and ([7|) we assumed that the electric field E in the BLG region was exactly 
the external one, E^xt- There is, however, an obvious additional contribution: the external 
electric field polarizes the BLG, inducing some charge asymmetry between the two graphene 
layers, which in turn give rise to an internal electric field, Eint, that screens the external one. 



V = {n/riQ — 2)e^nod/(2e,.eo)- 



(6) 



V = {2- n/n^)e^n^d/{2er£o). 



(7) 
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Figure 2. (a) y vs density for 
the epitaxial biased BLG device. 
Experimental data from Ref. |3] is 
shown as O . The line in the inset 
shows the unscreened result given 
by Eq. (b) Gap vs density for 
the exfoliated biased BLG device 
with chemical doping set to ng. 



To estimate Eint we can again apply a parallel plate capacitor model. The internal electric 
field due to the charge asymmetry between planes may thus be written as Ei^t = eAn/(2ereo), 
where — eAn is the induced charge imbalance between layers, which can be estimated through 
the weight of the wave functions in each layer, 

An = m - na = ^ i\^^A\,\.? + l<^Bi,kl^ " Wa2^ " \'fm}.?)^ (8) 

^ j,i=± k 

where the factor 2 comes from spin degeneracy, A^c is the number of unit cells and Aq = a^\/3/2 
is the unit cell area, jl is a band label, and the prime sum runs over all occupied k's. Note that in 
order to calculate An we must specify V, as it determines the amplitudes on the right hand side 
of Eq. dS]). On the other hand. An determines Eint, and in its turn Eint enters Eq. ^ through 
E to give V. Thus, a self-consistent procedure must be followed to calculate the screened energy 
difference between layers V. In particular, for the two experimental realizations of biased BLG 
discussed in Sec. 13. H the self-consistent equation that determines V reads 

V =[n-2no + An(n, V)]e^d/{2ereo) and V = [2n^ -n + An(n, V)]e^d/{2£reo), (9) 

respectively for exfoliated BLG |5j and epitaxial BLG [3] devices. It is worth noting that the 
screening correction expressed in Eq. ([9]) leads to a term in the Hamiltonian which is exactly 
the Hartree correction due to the charge asymmetry between layers. This approach has been 
followed in Refs. in the continuum approximation, and has been tested against ab initio 

calculations at half- filling in Ref. |16j . 

In Fig. [2{a) we compare Eq. ^ (right) with experimental results from Ref. [3] already 
mentioned in Sec. 13.11 Clearly, the self-consistent V given by Eq. ^ for e,. = 1 is a much 
better approximation than the unscreened result of Eq. ([7]) (see inset). In Fig. [SJb) we show the 
gap Ag as a function of carrier density n for the exfoliated biased BLG device, with realistic 
values of chemical doping no [5]. The gap is given by Eq. (jl]), with t± ~ 0.22 eV [5] and V 
obtained by solving self-consistently Eq. Q (left) for = 1. Note that for E^xt = we always 
have Eint = (the charge imbalance must be externally induced), and therefore we also have 
y = and Ag = 0. For this particular biased BLG device the present model predicts E^xt = 
for n = 2no, which explains the asymmetry for Ag vs n shown in Fig. [2jb). Note also that 
this device effectively provides a tunable gap semiconductor, as implied by Fig. [2^b): different 
gap values are achieved by tuning the gate voltage Vg, which controls the carrier density n. 
Note that the maximum A^ occurs when Vg reaches its maximum, which occurs just before the 
breakdown of SiOa. The breakdown field for SiOa is > IV/nm, meaning that Vg values as high 
as 300 V are possible for the exfoliated BLG device. As shown in Sec. 13. H Vg ~ ±300 V implies 
n — riQ^ ±22 x 10^^ cm^^, and therefore Fig. [2]|^b) nearly spans the interval of possible densities. 
It is apparent, specially for no = 5.4 x 10^^ cm~^, that when the maximum allowed densities are 
reached the gap seems to be approaching a saturation limit. This saturation is easily identified 
with the plateau shown in Fig. [^d) for Ag vs V , occurring for F > tj.. We may then conclude 
that this device enables the entire range of allowed gaps (up to tj_) to be accessed. 



4. Edge states in bilayer graphene 

The study of edge states in AS— stacked BLG given here is based on the ribbon geometry 
with zigzag edges shown in Fig. [3ja). Zigzag edges break translational invariance along their 
perpendicular direction, enabling us to write an effective one-dimensional Hamiltonian for a given 
momentum ka G [0, 27r[ along the ribbon. The first nearest-neighbor tight-binding Hamiltonian 
[Eq. ([1])] can then be written as 

2 

Hk = -tY,Yl 4;k,n(-^'^^''^^^kbi-k,n + 6i;fc,„-l) - t±Y^ 4;k,nhk,n + h.C, (10) 
i=l n n 

where aj;A:,n {bi-k,n) is the annihilation operator at momentum k and position n across the ribbon 
in sublattice Ai (Bi), i = 1,2 is the layer index and Dk = — 2 cos(A;a/2). 



4.I. Semi-infinite bilayer graphene 

Edge states in BLG are investigated by solving the Schrodinger equation, Hf^ 1?/^^) = E'^ IV'/c)- 
An analytic solution is possible for the semi-infinite BLG [17J, the same being true for SLG 
|12j . The wavefunction li/j^) is written as a linear combination of the site amplitudes along 
the edge's perpendicular direction, 1-0^.) = X^n X]i=i['^i(^' {o-i^kju) + Pi{k,n) \hi,k,n)], where 
we have introduced the one-particle states \ci,k,n) = c\.^^ |0), with q = ai,hi. Assuming all 
edge atoms belong to the A sublattice (without loss of generality), we require the boundary 
conditions ai{k,n — )• oo) = aj(/c,— 1) = f3i{k,n — oo) = /3i(/c,— 1) = 0, with i = 1,2, accounting 
for the existence of the edge at n = 0. Within our model, the Fermi energy of BLG always 
occurs at zero energy. Therefore, we expect zero energy edge states to have interesting physical 
consequences, and we set Ef^ = 0. As a result, the two sublattices become completely decoupled, 
and only the sublattice to which edge atoms belong can support edge states. This means that 
we always have /3i{k,n) = P2{k,n) = 0. 

With H}^ as in Eq. (jlOp . it is easy to write iV'fc) = as a transfer matrix that relates 
Qj(/c, n + 1) with ai{k, n), for i = 1,2. Then, applying the abovementioned boundary conditions 
we find that bilayer graphene supports two types of zero energy edge states localized at zigzag 
edges for 27r/3 < ka < 47r/3: one type restricted to a single layer and coined monolayer family, 
with amplitudes equivalent to edge states in SLG, 

ai{k,n) = Q and a2{k,n) = a2{k,^)Dle-'^'^; (11) 
and a new type coined bilayer family, with finite amplitudes over the two layers, 

Qi(A;,n) = ai(A;,0)L>^e-*^" and a2(A;, n) = -ai(A;, 0)Z)^-^ye-*^("-i) - , 

' (12) 

where the normalization constants are given by |a2(A;,0)p = 1 — and |ai(/c,0)p = 
(1 — DI)^/[{1 — + t\/t^]. An example of the charge density associated with Eq. ([12]) 
is shown in Fig. EKb) for = 0.2t, where the |ai(/c,n)p dependence can also be seen as the 
solution given by Eq. (jlip for \oi2{k,n)'\^ , apart from a normalization factor. 

We note that the same value of A, the penetration depth, is obtained from Eqs. (jlip and (jl2p : 
A = —1/ In \ Dk\. Nevertheless, the solution given by Eq. (jl2p has a linear dependence in n which 
enhances its penetration into the bulk. We expect these states to contribute more to self doping 
then the usual SLG edge states [13], as the induced Hartree potential which limits the charge 
transfer between the bulk and the edge will be weaker. 




Figure 3. (a) Bilayer graphene ribbon with zigzag edges, (b) Charge density for bilayer edge 
states at ka/2'K = 0.36 (top) and ka/lir = 0.34 (bottom), (c) Energy spectrum for a bilayer 
ribbon with zigzag edges, = 400, i_L = 0.2i: V = Q (top) and V = (bottom). 



^.2. Bilayer zigzag ribbons 

So far we studied locahzed states at the semi-infinite BLG. Experimentally, however, the relevant 
situation is a BLG ribbon. The band structure of a BLG ribbon with zigzag edges and 
unit cells in the perpendicular direction is shown in the top panel of Fig. [3]^c), obtained by 
numerically solving Eq. (jlOp . There are four partly flat bands at ii^ = for k in the range 
27r/3 < ka < 47r/3, corresponding to four edge states, two per edge. Strictly speaking, the 
edge states given by Eqs. ([TT]) and (fT2]) . and the other two localized at the opposite edge, are 
eigenstates of the semi-infinite system only. In the ribbon the overlapping of four edge states 
leads to a slight dispersion and non-degeneracy. 

The fact that edge states do exist in BLG ribbons imply several interesting physical properties. 
First we note that edge states induce a strong peak in the density of states (DOS) at zero 
energy, due to the fiatness of the energy bands. This DOS peak can be easily obtained as 
follows. Assuming we are sufficiently far from the Dirac points, and taking into account 
that t±/t <C 1 holds, we can write for the overlapping between monolayer and bilayer edge 
state families |T|'| ~ t\Dk\^, while the overlapping between two edge states of the bilayer 
family can be written as |T^^| ~ t±N iDkl'^^^ (the overlapping between two edge states of the 
monolayer family being identically zero). The energy dispersion may then be approximated by 
± \Dk\^ [(4*2 + t^N^/Diy/^ ± t^N/ \Dk\], which, near A;a ~ vr and assuming A^ > 1, can 
be further simplified as 

\±t^N\ka-T:f 
^'^\±t^\ka-^f /{t^N) ■ ^'^^ 

Thus the DOS induced by edge states has the form p{E) ~ E^'^^t^ /N , which is also found in 
SLG. Therefore, as predicted for SLG [18], we expect a Curie-like temperature dependence for 
the Pauli paramagnetism in BLG due to edge states - twice as large due to the presence of 
twice as many edge states in BLG. Furthermore, edge magnetism (close to zero temperature) 
was found to be possible in BLG due to edge states [19], similarly to SLG [12J. 

Of particular importance is the effect of a perpendicular electric field applied to the zigzag 
bilayer ribbon. The semi-infinite biased system has only one edge state given by Eq. (jlip . as 
the edge state having finite amplitudes at both layers [Eq. ()12p ] is no longer an eigenstate. In 
the bottom panel of Fig. [Sjc) we show the band structure of a BLG ribbon for V = t±/2. Two 
partially fiat bands for k in the range 27r/3 < ka < Att/2> are clearly seen at E = These 
are bands of edge states localized at opposite ribbon sides, with finite amplitudes on a single 
layer [Eq. (jlip ]. Also evident is the presence of two dispersive bands crossing the gap, which are 
reminiscent of the bilayer family of edge states. This dispersive states appearing inside the gap 



may contribute to the finite spectral weight recently observed using ARPES [3]. 
5. Conclusions 

We have studied the electronic behavior of BLG using the minimal tight-binding model that 
describes the system. Particular focus has been given to the presence of an external electric 
field perpendicular to the BLG system - biased bilayer - which gives rise to a finite gap in 
the spectrum, whose size is completely controlled by the applied voltage. The effect of the 
perpendicular electric field has been included through a parallel plate capacitor model, with 
screening correction at the Hartree level. The biased BLG thus realizes the first tunable gap 
semiconductor - a proof of principle regarding real applications of BLG. We have also addressed 
the effect of zigzag edges in BLG. We have found that BLG, as its SLG counterpart and other 
graphene based materials |20) . possesses zero energy surface states which can be divided into 
two families, giving rise to twice as many zero energy bands as in SLG. In the biased case half 
of the bands become dispersive inside the gap. 
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